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We study the (spin-)Peierls transition in quasi-one-dimensional dis- 
ordered systems, treating the lattice classically. The role of kinks, in- 
duced thermally and by disorder, is emphasized. For weak interchain 
interaction the kinks destroy the coherence between different chains at 
a temperature significantly lower than the mean-field Peierls transition 
temperature. We formulate the effective Ising model, which describes 
such a transition, investigate the doping dependence of the (spin-)Peierls 
transition temperature and discuss several implications of the picture de- 
veloped. The results are compared with the properties of the spin- Peierls 
system CuGe03. 
PACS numbers: 64.70Kb, 75.10Jm 



I. INTRODUCTION 

The tendency towards the lattice instability is strongly enhanced in quasi-one- 
dimensional systems. Conducting chains are inherently unstable against the periodic 
lattice distortion which opens a gap in the spectrum of electron excitations (Peierls 
instability) .0 TaS3, NbSe3, K0.3M0O3 are the examples of the quasi-one-dimensional con- 
ductors, which at a sufficiently low temperature become Peierls insulators.! The peri- 
odic distortion was also observed in organic (TTF-CuBDli, TTF-AuBDtS (MEM)- 
(TCNQ) 2 §, etc.) and inorganic (CuGe0 3 i, NaV 2 O 5 0&) compounds, consisting of weakly 
interacting spin chains. 

The excitation spectrum in the quasi-one-dimensional electronic and spin Peierls ma- 
terials is usually interpreted in terms of the excitations in isolated chains. Nonetheless, the 
(spin-) Peierls transition is a three-dimensional phenomenon, because electronic or spin 
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chains are immersed in a crystalline lattice of a bulk material. In the standard description 
of the (spin-)Peierls transition the lattice is treated in the mean-field approximation. 
In this approach, the excitations, which eventually destroy the ordered Peierls phase at 
a critical temperature, are the electron-hole (spin) excitations in isolated chains with a 
fixed lattice distortion. 

It is well-known, however, that in these systems there exist also excitations of a differ- 
ent nature: solitons, or kinks, which are the mobile domain walls. In this contribution we 
consider the Peierls systems, in which the lattice in the ordered phase is dimerized, i.e. 
long and short bonds alternate. In this case the distortion wave is commensurate with 
lattice and cannot slide as a whole. On the other hand, kinks, separating domains with 
the opposite sign of the dimerization, can propagate along chains. Solitons are mixed 
electron-lattice (or spin-lattice) excitations. Their topological nature is responsible for 
the anomalous charge-to-spin ratio, which was used to explain the unusual bevavior upon 
doping of the electronic Peierls material trans-polyacetylene. For spin-Peierls systems, 
the soliton picture is much less popular, although it was discussed also in this context in 
Ref. [Tl 

In addition to kinks, excited thermally or optically, there can be also kinks induced 
by disorder. Here we consider only one kind of disorder, which effectively cuts chains 
into finite segments. The interruption of chains may be caused by strong conformational 
disorder (e.g. bond bending in trans-polyacetylene) , or by substitution of atoms in chains 
(as in CuGeOs). The minimal energy lattice configuration strongly depends on whether 
the number of atoms in a chain is even or odd: an odd chain contains one kink. 00 while 
an even chain does not. The disorder- induced kinks are neutral and have spin |. In fact, 
any kind of the off-diagonal disorder results in the appearance of kinks,0'0 so the picture 
we consider here is rather general. 

Although kinks provide a coherent explanation of a vast variety of experimental data 
on optical and magnetic properties of trans-polyacetyleneJiMl their role in the Peierls 
phase transition was not properly studied. That may be related to the fact, that trans- 
polyacetylene has a large gap E g ~ 1.8eV, so that the transition into the undimerized 
state cannot be observed. On the other hand, the critical temperature in spin-Peierls 
materials is rather low (2-35 K), and there exist a very detailed information on these 
transitions. An important property of the recently discovered inorganic spin-Peierls com- 
pound CuGeOs is its ability to accept dopands, so that the disorder can be introduced 
into this system in a controlled way. This made possible the systematic study of the effect 
of disorder on the spin-Peierls transition.&"@ Finite chains in the spin-Peierls compound 
CuGe03 are formed upon the substitution of Cu atoms in the Cu02 spin chains by 
nonmagnetic Zn, which interrupts the spin-exchange. Similar effect is reached by the 
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substitution of Ge by Si.@ Doping by Zn or Si strongly decreases the Peierls transition 
temperature and gives rise to the antiferromagnetic phase,0H which coexists with the 
dimerization. 00 We will show, that these effects can be naturally explained in terms of 
the disorder- induced kinks. 

In this contribution we consider the (spin-) Peierls transition into a dimerized state in 
disordered quasi-one-dimensional systems, taking into account the soliton excitations. We 
believe that, due to their topological nature, kinks play an important role in the descrip- 
tion of this transition, even though the interaction between the chains can considerably 
modify their properties. The approach we propose is, in a sense, opposite to the mean- 
field description of the (spin-) Peierls transition: we neglect the electron(spin) excitations, 
as well as the small fluctuations of the order parameter, and leave only kinks. In other 
words, we assume that locally the dimerization persists well above the critical tempera- 
ture, and that the transition occurs when the coherence between the dimerization phases 
on different chains is established. The effective model, which describes such transition in 
clean systems, is the strongly anisotropic Ising model. We propose a simple modification 
of this model, which allows us to consider both thermally and disorder-induced kinks. 

Our approach is valid if the coherence of the (spin-) Peierls state is destroyed by the 
kinks at a temperature significantly lower than the mean field transition temperature, 
so that the gap in the spectrum of the (spin) electron excitations is still large at the 
phase transition point and all nontopological excitations can be neglected. As we shall 
see, this requires a rather weak interchain interaction. Another important approximation 
we make, is the classical treatment of the soliton motion. The classical treatment of 
lattice is justified, if the value of the relevant phonon frequency is much smaller than the 
value of the Peierls gap. This is usually the case for the electronic Peierls materials (e.g. 
trans-polyacetylene). On the other hand, in spin- Peierls materials with a relatively small 
spin gap this would require a particularly soft lattice, which, in principle, can happen 
due to a proximity to some structural transition.0 If, on the contrary, the Peierls gap is 
comparable to the optical phonon frequency, or is much smaller, then the consideration 
of the three-dimensional lattice has to be quantum. 



II. WEAKLY COUPLED FINITE CHAINS 

In this section we introduce a simple model of disordered quasi-one-dimensional (spin- 
) Peierls material. There are three essential assumptions on which our model is based: 
(i) interchain interaction is assumed to be sufficiently weak; (ii) off-diagonal disorder is 
assumed to break effectively chains into segments; (iii) lattice can be treated classically. 

Because of the weakness of the interchain interaction, in the first approximation chains 
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can be considered independently of each other. The typical length of a segment in 
a disordered system is inversely proportional to the concentration of impurities. The 
classical lattice configuration of a finite chain has to be determined by minimizing the 
total chain energy, which includes the electron energy (the energy of the spin system for 
spin-Peierls materials), as well as the lattice energy. Below we discuss the form of the 
minimal energy lattice configuration of finite chains. 

We begin by considering a conducting chain described by the Peierls-Hubbard (PH) 
Hamiltonian, 

H PH = - X] (*o + a(u n - u n+ i)) (c na c n+l(T + c ] n+la c ncT ) + U^2 c n ^c ni c nl 

n,cr n 

K 

+ -Un+i? , (1) 

^ n 

where the operator c na annihilates an electron with spin projection a at a site n, and u n 
is the shift of the ra-th atom in the chain direction. The energy of a half-filled chain with 
even number of atoms and periodic boundary conditions is minimal when the chain is 
uniformly dimerized, 

Un = (-) n U , (2) 

which corresponds to alternation of long and short bonds. The dimerization takes place 
at all values of U. For non- interacting electrons (U = 0),E!S 

oc e Xe p , (3) 

where X ep is the dimensionless electro n-pho non coupling constant, 

4a 2 

In the opposite limit of strong electron correlations, U 3> 4t , the value of the dimerization 
is determined by the interaction of the lattice with the low energy spin degrees of freedom, 
described by the spin-Peierls (SP) Hamiltonian, 

1 K 

H*P = 2 (^o + a\u n - u n+ i)) (S n ■ S n+ i) + — ^{u n - u n+1 ) 2 , (5) 



where 



In this case 
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Jo = — , and a = —a . (6) 



to V u 
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The number of bonds in even chain with periodic boundary conditions (a ring) is also 
even, so that there are two degenerate minimal energy lattice configurations, shifted by 
one lattice constant with respect to each other. The amplitude of the distortion wave, uq, 
for these two lattice configurations has opposite signs. The boundary conditions for finite 
chains produced by substitution of atoms are, however, free rather than periodic. For an 
open chain with even number of atoms there is only one minimal energy configuration, 
since the number of bonds in such chain is odd, so that the configuration, which has 
more short bonds, has lower energy (see Fig. la). The value of the dimerization changes 
somewhat from the middle to the ends of the chain, but the sign of the dimerization is 
everywhere the same. 

On the other hand, if the number of atoms in a chain is odd, the minimal energy 
lattice configuration has a form of a kink (see Fig. lb). The chain energy is practically 
independent of the position of the kink, provided that the distance from the kink to the 
nearest chain end is larger than the kink size. One unpaired spin is localized near the 
kink. Therefore, a spin | object in an isolated odd segment, contrary to a common belief, 
is localized not close to an impurity (chain end), but rather away from it. 

Next we consider the effects of interchain interaction. This interaction may origi- 
nate, for instance, from a relatively small hopping of electrons between chains (which 
corresponds to the exchange between spins sitting on neighbouring chains in spin-Peierls 
materials), long range Coulomb interaction between the electrons and deformations of the 
three-dimensional lattice. Independently of the origin, the interchain interaction tends 
to create a coherence between the phases of the distortion waves on different chains. 
In a clean material it determines a type of crystalline structure. For simplicity we will 
assume, that the lattice configuration has the lowest energy when the sign of the dimer- 
ization is the same on all chains. (Actually, one can show that in CuGe0 3 the distortions 
in neighbouring chains along 6-axis have opposite phases, but it does not modify our 
conclusions.) 

When a finite chain, created by disorder, is surrounded by other chains, its lattice 
configuration is different from that of an isolated chain. For instance, there appear two 
kinds of even chains: those in which the sign of the dimerization, favoured by the bound- 
ary conditions, is the same as the sign of the average dimerization in the bulk material, 
and those in which it is the opposite (see Fig. 2a). In the latter case the energy loss 
due to interchain interaction is proportional to the chain length. Therefore, despite the 
weakness of interchain interaction it will be energetically more favourable in a sufficiently 
long chain to create a kink near one end of the chain and antikink near the opposite end, 
so that the sign of the dimerization between the domain walls is the same as the sign of 
the averaged dimerization (see Fig.2b). Similarly, the interchain interaction binds kink 
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in an odd segment to one of the chain ends (see Fig.3). This, in a sense, restores the 
conventional picture of a spin ~ sitting relatively close to impurity at an average distance 
of the order of a kink size £ (see also Ref. 



III. QUALITATIVE DISCUSSION OF THE PHASE TRANSITION IN CLEAN 

AND DISORDERED SYSTEMS 

In this section we discuss qualitatively the (spin-) Peierls phase transition in quasi- 
one-dimensional systems with very weak interchain interaction. First we consider the 
clean systems. The long range order in an infinitely long isolated chain is destroyed by 
thermally excited kinks at any non-zero temperature. This happens because the energy of 
kink does not depend on its position in an isolated chain. In the presence of the interchain 
interaction the string with the opposite sign of the dimerization is created between the 
soliton and antisoliton. It costs energy proportional to the length of the string I, which 
results in a linear potential between kinks, 

V(l) = XI . (8) 

The coefficient A measures the strength of the interchain interaction and may be called 
the string tension. The potential tends to bind solitons and antisolitons into pairs. The 
phase transition now occurs at a finite temperature, when the thermally excited soliton- 
antisoliton pairs dissociate. This happens when the pair size, 

R(T) = | , (9) 

becomes comparable to the distance between the thermally induced kinks in an isolated 
chain, 

d(T) = exp(|) (10) 

where fi is the kink excitation energy, which is of the order of gap in electron (spin) 
excitation spectrum. In this way for weak interchain coupling, A <C fi, we obtain approx- 
imately 



T c (0) ~ ^ • (11) 



t 

A 

In the doped Peierls system, apart from the thermally induced kinks, there exist also 
kinks induced by disorder, whose density is proportional to the concentration of dopands. 
The total density of kinks is 
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Iltot = ntherm + n dis = 6 t + Cx , (12) 

where C is some numerical factor of the order of 1. We can now estimate the phase 
transition temperature T c (x) in the same way as we did for a clean system, using 1/ritot 
as the average distance between the kinks. At small x, when the density of the disorder- 
induced kinks ridis is much smaller than the density of the thermally excited kinks n t h erm , 
the phase transition temperature decreases linearly with x, 

T c {x)=T c {0)(l-Ax) , (13) 

and the coefficient A is large, 

- " (14) 



A more careful estimate of A will be given in the next section. If, on the other hand, 
ndis n t herm, the phase transition temperature becomes inversely proportional to the 
concentration, 

T c (x) = - . (15) 

Ob 

The interpretation of the last equation is, that at the phase transition point the typical 
distance from kink to the nearest end of odd segment, ~ T/X, becomes comparable with 
the average length of the segment ~ 1/x. This result is not applicable at very large x: 
when an average length of the segment becomes of the order of soliton size, 1/x ~ £ i 
there will be no ordering even at T = 0. 

Because the density of the thermally induced kinks is exponentially small at T <C /i, 
the disorder-induced kinks begin to dominate at a very low doping concentration, 

*c~-1h(y) » ( 16 ) 

and the crossover between linear and | decrease of the transition temperature is very 
sharp for A <^ ji. 



IV. EFFECTIVE ISING MODEL 

The model qualitatively discussed above, which describes the statistical properties 
of kinks, induced thermally or by disorder in the ensemble of weakly interacting finite 
chains, can be formulated as an Ising-type model. To this end we divide chains into 
cells containing two bonds. The sequence of bonds in each cell can be either short-long 
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or long-short (positive or negative dimerization), which is decribed by the Ising variable 
a — ±1. The uniformly dimerized lattice corresponds to the ferromagnetic ordering 
of the Ising spins. Kinks are domain walls in Ising variables, i.e., o n = —1 for n less 
than some m and a n = +1 for n > m, or vice versa. In this model one neglects the 
actual size of the kink, which is reasonable as long as it is much smaller than the average 
distance between the kinks. The weak interaction between different chains we treat in 
the mean-field approximation.il The energy of the chain containing N Ising spins has the 
form, 

a N N a' 

E N =^E( 1_ °n°n+l) - h °n - ~~ (^1 + + fi' . (17) 

Z n=l n=l Z 

Here the first term describes the energy cost of a kink (equal to fi) , while the second term 
describes the interaction with the mean field h of neighbouring chains. The third term 
is introduced to model different kinds of finite chains in (spin-)Peierls material. It has a 
form of the interaction with the "magnetic field" Hl = ufj,', applied at the left end, and 
Hr = u'fi' at the right end. Here v and v' take values ±1, which corresponds to four kinds 
of finite (spin-)Peierls chains: even segments in phase {y — v' — +1) and out of phase 
{y = v' — —1) with an average dimerization, and odd segments with one end in phase 
and other out of phase [y = +1, v' = —1 or vice versa). In the latter two cases, the 
opposite orientation of the magnetic field at the chain ends creates a kink in the lowest 
energy configuration of the Ising chain, provided that // > fi. In what follows, we will 
restrict ourselves to the case of infinite jj . 

The partition function of the Ising chain can be easily found by means of the transfer 
matrix method.0 For v' = v the result is, 



N 



e -uh 



2 

while for v' = —v one has, 



( A iV+i + \ N+1 ) + z/sinV(A^ +1 - A^ 1 )! , (18) 



Z%- = V = ^(A^ 1 - A^ 1 ) . (19) 



Here we denoted by A± the eigenvalues of the transfer matrix, 



A± = cosh h±J (sinh h) 2 + e~ 2 ^ , (20) 
and if) is defined by, 

sinh h , s 

sm^= ; _ = (21) 

(sinh h\ + e~ 2lM 
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Above we used the notations: h = (3h and p, = (3p. 

The averaging over disorder reduces in our model to the averaging over N, u, and v' . 
For simplicity we assume that for each N all the four types of finite chains are possible. 
The distribution of the chain lengths Pn for small concentration x of randomly positioned 
dopands (e.g. Zn in CuGeOs) is 

p ~ = M exp h£>) ' (22) 

where the average chain length is related to the concentration of dopands by, 

(N) = -, (23) 

x 

(factor 2 appears because the length of the unit cell in the effective Ising model is equal 
to 2a). For small concentration x one can substitute the summation over N by the 
integration. 

The averaged free energy per unit cell is 

</> = tAa E P * F n , (24) 

^\ iV / N,v,v> 

where Fff = —T In Zpf is the free energy of a finite chain. The mean-field value h is 
proportional to the averaged value of the Ising spin (the average dimerization), 

h = X(a) = -A^ , (25) 

where A is the coupling constant. 

The transition temperature, T = T c , is obtained by soving equation!! 

d 2 (f) 

-^r(^ = o) = i, (26) 



Evaluation of the second derivative of the average free energy gives 
A ' " [ ^C(2, J) - 1 - 26 7 + 4e 7 3 ( /3( 7 ) - i- 



where ((z, q) is the Riemann zeta function 



27, 



(27) 



and f3(q) is,§ 



OO 1 

C( " 9) =S^- (28) 



«')=S^- (29) 
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We also introduced the dimensionless parameter 

7 = xe^ , (30) 

which is, roughly speaking, the ratio of the densities of the disorder-induced and the 
thermally-induced kinks, and 

e = -e* In fcoth . (31) 



2, 

For T <C n e is close to 1. 

If now 7< 1, ie. the average distance between the thermally induced kinks is much 
smaller than the average chain size, the equation for T c (x) can be approximately written 
as follows: 



w(«fo) U_ 2(XW (jl 



r. ( «) v- 2exexp m)) = u (32) 

At x = the critical temperature is determined by the equation, 

exp lm) = — • (33) 

which coincides exactly with the condition R(T c (0)) = d(T c (0)), which we used in the 
previous Section (see Eqs. (0), fllCf), and ([TT]) ) . The magnitude of the linear slope of 
T c (x)-curve at x = (cf. Eqs.(0) and QI|])) is 

T c (0) Aln(^ X ' 
In the opposite limit 7 ^> 1, when the disorder-induced solitons dominate, we obtain 



r e( «) = I , (35) 



from Eq.(|27|), 



which has to be compared with Eq. (|T5|) of the previous Section. 

Eq.(|27|) was solved numerically for the concentration of dopands x < 0.08 and several 
values A. The results are presented in Fig. 4. The value of \i for each A was chosen to 
maintain T c (0) = 14X. We see that the crossover between the linear and ^ behavior of 
T c (x) is indeed very sharp. 



V. CONCLUSIONS 

The treatment carried above shows the important role played by solitons in describing 
properties of both clean and doped (spin-)Peierls systems. For the electronic Peierls ma- 
terial like polyacetilene, there are already a lot of experimental evidences of the existence 
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of kinks (although due to a large kink energy they do not contribute significantly to the 
thermodynamic properties of this material at room temperature). 

The mean-field soliton lattice solution was used to describe the thermodynamics of 
the incommensurate phase of spin-Peierls systems in strong magnetic field.Q0Hil We 
have shown, however, that solitons could strongly affect the properties of spin-Peierls 
systems, both clean and doped, also at H = 0. The small parameter \ j \x (where A is the 
strength of the interchain interaction and /i is the kink energy of the order of spin gap) is 
responsible for the sharp suppression of T c by doping (cf. Eqs. ffLU) and ([34]) ). Qualitative 
estimates show, that the SP phase is completely suppressed, when the average distance 
between impurities becomes of the order of the coherence length of the SP ordering £, 

I ~ £a , (36) 

or, in other words, when the average segment length becomes of the order of the kink 
size. For CuGeC>3, £ ~ 8-10 lattice spacings, so that the critical concentration for 
the complete suppression of SP transition x c cannot, in any case, exceed ~ 10-12 %. 
Experiments give a somewhat smaller value: x c ~ 7-8 % for Zn and x c ~ 2-3 % for Si. 
The extra "efficiency" of Si, as compared to Zn, is, most probably, related to the fact 
that Si, substituting Ge, is located between two CuCVchains and thus influences two 
chains simultaneously. Therefore, two kinks are created per one Si and only one kink 
per Zn. Correspondingly, Si doping may be expected to be two times more efficient in 
suppressing T c than the doping by Zn. 

An additional argument in favour of solitons induced in doped (spin-)Peierls systems, 
comes from certain analogy between the properties of doped CuGeO"3 and the behaviour 
of this material at strong magnetic fielcS (in which case the soliton picture was checked 
experiment allyEHJHI) . The notion of solitons can, therefore, provide a unified approach to 
this problems. 

Yet another problem, appearing in the discussion of the doped spin-Peierls systems, 
is the interplay between the SP and antiferromagnetic ordering. As was first observed in 
Refs. ^1] and [24], the two phases coexist in some interval of the concentration of dopands. 
Soliton picture provides a natural explanation of this phenomenoni! As was discussed 
above, the kinks created by doping carry spin |. They induce in their vicinitiy the 
antiferromagnetic correlations, which decay at the length scale £ along the chain (see 
Fig. 5). The antiferromagnetic microregions, created in this way, interact with each other, 
giving finally the long-range magnetic order, but with a reduced sublattice magnetization, 
spatially varying due to disorder. It is important that the spins of the separated kinks 
are free, so that they can ajust to the mean field created by other spins. In other words, 
the disorder, inducing unpaired spins, does not introduce frustration in the magnetic 
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interaction. Once the coherence between the spins is established, the Neel temperature, 
T^, should grow, roughly speaking, linearly with the density of the disorder-induced 
spins, i.e., Tjy oc x at small x. The detailed treatment of the interplay of the SP and 
antiferromagnetic ordering will be given in a separate publication. 

A similar picture, but without an apparent use of the soliton concept, in which the 
local antiferromagnetic correlations are created in the regions where the lattice ajusts to 
impurities, was suggested in Ref. Creation of free spins by nonmagnetic impurities in 



the rigidly dimerized chains was considered in Ref. |36 



Our model is, of course, oversimplified, as we assumed that the phase transition 
into disordered state occurs entirely due to the loss of coherence between the phases of 
the order parameter on different chains. In practice, the conventional spin and phonon 
excitations have to be included to get a reasonable description of the excitation spectrum, 
magnetic susceptibility and other properties of spin-Peierls materials. The aim of this 
contribution, is, however, to show that the large fluctuations of the order parameter 
(kinks), which are neglected in the mean-filed approach, can be important in the quasi- 
one-dimensional systems. 

Summarizing, we can say that the concept of solitons, which is well established for 
the electronic Peierls materials, is also very useful for the discussion of spin-Peierls sys- 
tems. It gives a natural framework for the description of the properties of disordered 
Peierls systems and allows one to explain the dependence of the spin-Peierls transition 
temperature on doping, the interplay of the spin-Peierls and antiferromagnetic ordering, 
the simultaneous influence of doping and magnetic field, etc. 
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FIGURE CAPTIONS 



FIG.l. Schematic view of the minimal energy lattice configuration of an isolated 
finite (spin-)Peierls chain with even (a) and odd (b) number of atoms (double and single 
lines correspond, respectively, to short and long bonds). The order parameter (thin line) 
changes sign in odd chain, and one unpaired spin (indicated by an arrow) is located near 
the kink. 

FIG. 2. Finite segment (thick lines) with even number of sites out of phase with the 
surrounding chains (thin lines) creates a string of unfavourable phase (a). The string can 
be removed by creating a kink and antikink near the chain ends (b). 

FIG. 3. Kink in an odd segment moved to one of the chain ends to minimize the 
interchain interaction energy. 

FIG. 4. Dependence of the (spin-)Peierls transition temperature on concentration of 
dopands for A = 3,2, 1,0.5, O.IK (the curves, respectively, (a) - (e)). 

FIG. 5. Antiferromagnetic correlations induced in the vicinity of kink. 
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